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Abstract 



Gravitational waves generated by inspiralling compact binaries are investi- 
gated to the second-post-Newtonian (2PN) approximation of general relativ- 
ity. Using a recently developed 2PN-accurate wave generation formalism, we 
compute the gravitational waveform and associated energy loss rate from a 
binary system of point-masses moving on a quasi-circular orbit. The cru- 
cial new input is our computation of the 2PN-accurate "source" quadrupole 
moment of the binary. Tails in both the waveform and energy loss rate at 
infinity are explicitly computed. Gravitational radiation reaction effects on 
the orbital frequency and phase of the binary are deduced from the energy 
loss. In the limiting case of a very small mass ratio between the two bodies 
we recover the results obtained by black hole perturbation methods. We find 
that finite mass ratio effects are very significant as they increase the 2PN 
contribution to the phase by up to 52%. The results of this paper should 
be of use when deciphering the signals observed by the future LIGO/VIRGO 
network of gravitational- wave detectors. 
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I. INTRODUCTION 



The LIGO/ VIRGO network of kilometer-size interferometric detectors of gravitational 
waves is expected to be in operation by the turn of the century (see for a recent 
review). The most promising targets for this network are the gravitational waves emitted 
during the radiation-reaction-driven inspiral of binary systems of compact objects (neutron 
stars or black holes). Crucial to the successful detection and deciphering of such waves will 
be the availability of accurate theoretical templates for the inspiral gravitational waveforms 
0-^. Much theoretical effort is currently being spent on developing improved formalisms 
tackling the generation of gravitational waves by general material sources and/or on applying 
existing formalisms to the explicit computation of increasingly accurate inspiral waveforms. 
Among the existing generation formalisms, the one proposed by us can, in principle. 



be developed to an arbitrarily high accuracy. Recent work by one of us [|I0| has succeeded 
in pushing its accuracy to the second post-Newtonian (2PN) level, i.e. in deriving general 
expressions for the asymptotic gravitational waveform, as a functional of the matter distri- 
bution in the source, which take into account all contributions of fractional order beyond 
the leading ("quadrupole formula") term. Here e ~ v/c ~ [Gm/rc^y/'^ denotes the small 
parameter entering the post-Newtonian expansion appropriate to the description of slowly 
moving, weakly stressed, weakly self-gravitating systems. 

The object of the present paper is to apply the 2PN-accurate generation formalism of 
Refs. [0-0] to the specific case of an inspiralling compact binary. This is a non trivial 
task as the end results of Ref. |jTO| contain some complicated three-dimensional integrals 
which are mathematically defined by a procedure of analytic continuation. Note that this 
contrasts with the end results of our previous 1.5PN-accurate generation formalism f^HH 
which contained only integrals extending over the (compact) support of the material source. 
The higher complexity of the 2PN-level is due to the appearance of terms associated with the 
cubic nonlinearities of Einstein's equations. We shall show below how to explicitly compute 
these terms in the case of binary systems. 

Theoretical waveforms such as the one computed below are useful to define parametrized 
"chirp" templates to be cross correlated with the outputs of the LIGO/ VIRGO interfero- 
metric detectors. For this technique to be successful, the templates must remain in phase 
with the exact general-relativistic waveform as long as possible. The phase of the signal is 
determined by the rate of change of the orbital period resulting from gravitational radiation 
reaction effects. Following the usual heuristic approach (which has been validated in detail 
at the leading order |11|), the effect of gravitational radiation reaction on the orbital period 
can be computed from the losses of energy (and angular momentum) |jl2|] at infinity. We 
shall therefore pay special attention to the computation of the 2PN-accurate energy loss 
from a (quasi-) circular compact binary which is obtained here for the first time, together 
with the resulting orbital phasing of the binary. 

Several investigations have recently focused on the computation of the energy loss and 
waveform in the case of a very small mass ratio between the two bodies [|T^-|TB[. Notably 
the energy loss has been computed in this case both numerically |T5[ and analytically 
up to an order going well beyond the 2PN-level. The test-mass limit of our result agrees 



with the 2PN truncation of the results of Refs. |]T5|,|T6[. However, we find that finite mass 



effects change very significantly the 2PN numerical contribution to the accumulated orbital 
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phase. Our main results, completed by the contributions due to the spins of the orbiting 



bodies, have been briefly reported in [T^ • Note that our formula for the energy loss has been 
conflrmed by an independent derivation based on a different, albeit less rigorous, method 
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The organization of the paper is the following. In Section |I| we write down the results 
of the 2PN generation formalism in a form convenient for the application to inspiralling 
binaries and we outline our strategy. In Section |IT^ we compute the 2PN-accurate "source" 
moments of mass-type (especially the quadrupole moment £ = 2) in the case of a binary 
system made of two point masses. In particular a crucial cubically nonlinear term is obtained 
in the subsection [111 C| . The expressions for all the relevant source moments are given in 
Section which deals with the explicit computation of the 2PN-accurate waveform and 
energy loss rate (including relevant tails). The instantaneous orbital phase of the binary is 
computed at the end of Section Technical details are relegated to several appendices: the 
conserved mass monopole and dipole moments are considered in Appendix Appendix ^ 
presents an alternative derivation of the cubically nonlinear contribution to the quadrupole 
moment, which is valid for N-body systems; and Appendix y is a compendium of various 
formulas for moments. 



II. SUMMARY OF THE 2PN-ACCURATE GENERATION FORMALISM 

Let us flrst recall that in a suitable "radiative" coordinate system X^^ = {cT, X*) the 
metric coefficients, say Ga/3{X''), describing the gravitational field outside an isolated system 
admit an asymptotic expansion in powers of R~^, when R = |X| —>■ oo with T — R/c and 
N = X/i? being fixed ("future null infinity"). The transverse-traceless (TT) projection 
of the deviation of Gais^X'^) from the fiat metric (signature — 1,+1,+1,+1) defines the 
asymptotic waveform = {Gkm{^) — ^km)'^'^ (latin indices i, j, k,m. . . range from 1 to 3). 
The 1/R part of can be uniquely decomposed into multipoles: 



AG °° 1 r 

~j^^N^L-2ea,i^j)tL^2{Tn)} + O (^) • (2.1) 



The "radiative" multipole moments Ul and Vl (defined for £ >2) denote some functions of 
the retarded time Tr = T — i?/c taking values in the set of symmetric trace-free (STF) three- 
dimensional cartesian tensors of order i. Here L = ii . . .ie denotes a spatial multi- index of 
order £, Nl-2 = Ni^ ■ ■ . Ni^_^, X(^ij) = \{X^ + X^i) and 

^.....(N) = (fa - iV,iV.)(i,„. - iV,iV„) - l(fa - iV.iV,)(4„. - iV.iV„) . (2.2) 



[For the convenience of the reader we summarize our notation in Ref. |T9|.] As indicated 
in Eq. (|2.1j ), for slowly moving systems the multipole order is correlated with the post- 
Newtonian order. The coefficients in Eq. ( p^.l[ ) have been chosen so that the moments Vj^ 
and reduce, in the non-relativistic limit c +oo (or e ^ 0), to the £-th time-derivatives 
of the usual Newtonian mass-type and current-type moments of the source. At the 2PN 
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approximation, i.e. when retaining all terms of fractional order with respect to the 

leading (Newtonian quadrupole) result, the waveform (|2.1| ) reads 
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— NabUijab + -£ab{iVj)acNbc 

1 2 ■ 

■^NabcUijabc + Y^£ab(iyj)acdNbcd 

^abcdUijabcd ~l~ gg ^ ab{i^j)acde-^bcde 



(2.3) 



The rate of decrease of the Bondi energy Eb with respect to the retarded time Tr = T — R/c 
is related to the waveform by 



dE 



B 



32nG 



dT, 



R 



R^dn{N) . 



(2.4) 



At the 2PN approximation this yields (with f/(") = d''U/dT]^ 
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A 2PN-accurate gravitational wave generation formalism is a method allowing one to 
compute the radiative moments entering Eqs. 



r^D and (2^) in terms of the source variables 
with an accuracy sufficient for obtaining the waveform with fractional accuracy 1/c^. The 
latter requirement implies, in view of Eq. (|2.3|) , that one should (at a minimum) compute: the 
mass- type quadrupole radiative moment Ui-^i^ with 1/c^ accuracy; the mass-type radiative 
octupole Ui-^i^i^ and the current- type radiative quadrupole Vi^i^ with 1/c^ accuracy; Ui^i^i^i^ 
and Vi^i^i^ with 1/c^ accuracy; Ui^i^i^i^i^ and Vi^i^t^i^ with 1/c accuracy; and Ui^i^i^i^i^i^ and 
V^jjjjgj^jg with the Newtonian accuracy. Note that these requirements are relaxed if one is 
only interested in getting the energy loss rate with 2PN-accuracy. In that case, Eq. ( p.5|) 
shows that one still needs Ui^i^ with 1/c^ accuracy, but that it is enough to compute Ui^i^i^ 
and Vijjj with 1/c^ accuracy, and Ui^i^i^i^ and Vi^i^i^ with Newtonian accuracy. 

In our generation formalism, the link between the radiative multipoles Ul and Vl and 
the dynamical state of the material source is obtained in several steps involving as inter- 
mediate object a certain vacuum "canonical" metric g^i^y{x^^-^) expressed in terms of some 
"canonical" multipoles and Sl (alternatively referred to as algorithmic moments in |P). 
On the one hand, the matching of S'l^t^l^can) a (PN-expanded) near- zone solution of the 
inhomogeneous Einstein equations allows one to compute and Sl in terms of some suit- 
ably defined "source" multipoles [source], J/, [source]. On the other hand, the computation 
of nonlinear effects in the wave zone allows one to compute Ul and Vl as functionals of Ml 
and 5*^. The final result for the 2PN- accurate generation formalism reads (when working 
in an initially mass-centred coordinate system, i.e. such that the canonical mass dipole Mj 
vanishes for all times) 
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for the moments that need to be known beyond the IPN accuracy, and 

UL{Tn)=li^\Tn)+0{e') , 



R) 
Vl{Tr) 



R) 
{Tr) 



+ 0{e' 



(2.6a) 
(2.6b) 
(2.6c) 



(2.7a) 
(2.7b) 



for the other ones. Eqs. ( p.6|) involve some integrals which are associated with tails; these 
integrals have in front of them the total mass-energy m of the source, and contain a quantity h 
which is an arbitrary constant (with the dimension of time) parametrizing a certain freedom 
in the construction of the radiative coordinate system (T, X) . More precisely, the link 
between the (Bondi-type) radiative coordinates = [cT, X*) and the (harmonic) canonical 
coordinates xff. 
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+ 0(£^) + 0(l/r, 
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(2.8) 



Except for the computation of Uij which requires the knowledge of the mass quadrupole 
source moment Jjj with 2PN accuracy, the computation of the other multipole contribu- 
tions to the waveform can be obtained from IPN-accurate expressions of the mass-type and 
current-type source moments which have been obtained for all values of i in Refs. and 

respectively, as explicit integrals extending only on the compact support of the material 
source. [Note that there are no 1/c^ contributions in the source moments.] Let us illustrate 
the structure of the IPN results by quoting the simpler IPN mass-type source moments: 



hit) 



^'"^^(''") + 2(2£ + 3)c2 dt^ 
4{2e+l)x,L 9(Ti(t,x" 



(£+l)(2£ + 3)c2 



dt 



+ 0(6' 



(2.9) 



The (compact support) matter densities appearing in Eq. ( |2.9| ), and their generalizations 
discussed below, are defined from the contravariant components (in the harmonic, "source" 
coordinate system x'^) of the material stress-energy tensor T'^'^ as 



cr(t,x) 

(Tij(t,X.) 



roo(t,x) + r^^(t,x) 



T°^(t,x) 

r^(t,x) 



(2.10a) 

(2.10b) 
(2.10c) 



The powers of c introduced in Eqs. ( |2.10D are such that the a's have a finite non zero limit 
as 1/c ^ 0. One associates to the matter densities ( p.lO|) some Newtonian-like potentials, 

say 
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[/(x,t)=G/^^a(x',t), 
J |x — x'l 

X(x, t)=G J rfV|x - xV(x', t) 
P,,(x,t) =G J 



dh 



X — X 



d,Ud,U - -6,,dkUdkU 



(x',t) 



(2.11a) 
(2.11b) 
(2.11c) 

(2.11d) 



[Note the equation satisfied by the X potentiah AX = 2U.] The 2PN-accurate source 
moment /l(^) has been obtained in Ref. |TD| and expressed in terms of the matter densities 
a, ai, aij, the potentials U, Ui, Pij and the trace P = P^s- The result (see Eq. (4.21) of 
Ref. [|TU1) is 



a + -{auU ~ aP) 



2c2(2£ + 3) 



c2(£+l)(2£ + 3) 
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ttGc^ V 4 



I 14 " 
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2(2£ + 1) 
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2(2£ + l)|xpx,L 
c4(£ + l)(2£ + 3)(2£ + 5) 



TXGC^ 



XL 



-P^jd,,U - 2UAdiU + 2d,UjdjUi - -{dtUf - Ud^U 



+ 0(£^) . (2.12) 



This expression obviously reduces to Eq. (|2.9|) at the IPN-level. It was also shown in 
Eq. (4.13) of Ref. [|IU| that the current moment Jl, which needs only IPN accuracy, can be 
written as 



|X| 3^L— l>a o2 

2c2(2£ + 3) 



{2i+l)xL-l>ac ^ 

'c^{i + 2){2£ + 3) * 



dkU{dbUk - dM) + -dtUdbU 



CTbc + 



AtiG 



dbUdJJ 



] + 0{e'). (2.13) 



This form is equivalent to the result previously derived in Ref. (Eq. (5.18) there). See 
below for comments on the symbol FPb=o in front of Eqs. ( 2.12|) - (|2.13|) . 

In this paper, it will be convenient to split the potential Pij of Eq. ( |2.11d[ ) into a compact- 
source potential Uij and a nonlinear potential Wij according to 



Pij = Uij - Wij + -SijWss , 



(2.14a) 



where Uij and Wij are defined by 



6 



[/,,(x,t) = G J 



-aij{^',t) , 

-[d,Ud,U]{^',t) . 



X — X' 



(2.14b) 
(2.14c) 



[We have AWij = diUdjU ] The compact-source potential Uij should not be confused with 
the previously defined radiative quadrupole moment Uij. Also note that the notation Wij is 
the same as used in Ref. ||10| to denote a different potential (which is a retarded version of 
the potential Pij) but which will not be used in this paper. Let us note for future reference 
that the trace of Wij as defined here satisfies 



$(x,t) 
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G 



[af/](x',t) 



(2.15a) 
(2.15b) 



The potentials introduced above are connected by the approximate differential identities 

dtU + diUi = 0{e^) , (2.16a) 

We replace in Eq. ( ^.12| ) the potential P^j by Eqs. ( ^.14| ) above, and dtU by the spatial 
derivative —diUi. It is also convenient to perform an integration by parts, using dkUd^Ui = 
^[A{UUi) — UAUi — UiAU], which is justified by a reasoning similar to the ones followed in 
Sect. IV of Ref. [0]. This leads to the form we shall use as starting point in this work: 
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d,Ud,U, - -d.Ud.U, 



dfai 



\dl{U') + W,,d,,U 



(2.17) 



For simplicity, we henceforth drop most of the post-Newtonian error terms as they are 
usually evident from the context. 

The symbol FP^^q iii Eq. ( p.l7|) stands for "Finite Part aX B = 0" and denotes a 
mathematically well-defined operation of analytic continuation. Let us recall its precise 
meaning (see [^,^| for details and proofs of the applicability of such a definition): one 
considers separately two functions of one complex variable B defined by the integrals: 
h{B) = /^^d3x|x|S/(x), I2{B) = Jy^d^:>c\x\^f{x.), where K is, say, the ball < |x| < tq 
and V2 the complementary domain: |x| > Vq. If the function /(x) is, say, continuous in 
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and has, at most, a polynomial growth 0(|x|^) when |x| oo, the integral defining Ii{B) is 
convergent if the real part of B is large enough, say Re{B) > —3, while the integral defining 
l2{B) is convergent if Re{B) < —p — 3. If moreover the function /(x) admits an asymptotic 
expansion ("multipolar expansion") of the form /(x) ~ Efc<p/fc j;^|x|'^n'^ as |x| — >■ oo (we do 
not include logarithms of |x| to simplify), the function I2{B) can be analytically continued 
as a meromorphic function in the complex B plane up to arbitrarily large values of Re{B) 
(with possible poles on the real axis). Finally, one considers the sum of Ii{B) and of the 
analytic continuation of l2{B) for values of B near 0: the constant term (zeroth power of 
B) in the Laurent expansion of Ii{B) + Analytic Continuation(/2(-B)) around B = defines 
FPb=o / c?^x|x|'^/(x). It is easily shown that this definition is independent of the choice 
of the intermediate radius tq used to split in two regions. Note that, in principle, one 
should introduce a length scale to a-dimensionalize |x| before taking its Bth power in ( p. 171 ). 



This is superfiuous in our case because []10| has shown that there arise no poles at 5 = 
and therefore no associated logarithms. 

Although the result ( p.l7|) is mathematically well-defined (contrary to the results of 



Refs. [p2| , p3| which are expressed in terms of undefined, divergent integrals), it is a non 
trivial task to compute it explicitly in terms of the source variables only. This will be done 
in the next section, in the case where the source is a binary system of non rotating compact 
objects (neutron stars or black holes). To this end, we shall heuristically represent the stress- 
energy tensor of the material source as a sum of Dirac delta functions 6. More generally, the 
stress-energy tensor of a system of (non rotating) compact bodies is formally given by 

r-(x,*)^E".ff^|^(x-y.W), (2.18) 

where rriA denotes the (constant) Schwarzschild mass of the Ath compact body. This yields 
for the source variables (|2.10| ) 

a(x, t) = Y: t^A{t) 1 + ^ U(x - y^(t)) , (2.19a) 

A=l \ / 

N 

ai(x, t) = J2 f^A{tyj{^ - yA{t)) , (2.19b) 

A=l 
N 

a,,(x,t) = fiAit)vAv'JiK - YAit)) , (2.19c) 



A=l 



where = dy\/dt and 

/iA(t) = mA[l + {d2)A + {dA)A] , (2.20a) 
d^ = -("^-v\ , (2.20b) 



d4 = ^ (^v^ + hj^^ - AU,v, - 2$ + + 4?7,,) , (2.20c) 
the notation V being a shorthand for the combination 

V^U+^d^X, (2.21) 
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which is the potential appearing naturally in the IPN near-zone metric in harmonic co- 
ordinates. The subscript A appearing in Eq. ( p. 20a ) indicates that one must replace the 
field point x by the position of the Ath mass point, while discarding all the ill-defined 
(formally infinite) terms arising in the limit x y^. For instance, 

iUU = G E ^"(;^(^ + "^/^'^ . (2.22) 



[Note that the second time derivative appearing in V, Eq. (p.21|) , must be explicated before 



making the replacement x — y^(t).] Although we do not claim to have verified it by a 
detailed proof, we feel secure that the formal use of delta functions can be justified at the 
2PN accuracy (and even at the 2.5PN accuracy) by combining our generation formalism 
with the results of Ref . . Indeed, the latter reference showed (by a matching technique) 
that the metric generated by a system of well-separated strongly-self-gravitating bodies was 
equal, up to the 3PN level, to the metric generated by a mathematically well-defined version 
of delta functions. See Section III.C of Ref. for a discussion, at the IPN level, of how to 
combine the two formalisms. 



III. THE 2PN-ACCURATE MASS MOMENTS OF A COMPACT BINARY 

It is convenient to split the starting formula ( p.l7| ) into three types of contributions, say 

/.=/r+/r+/r. (3.1) 

Here, ("Compact terms") denotes the terms where, because of the explicit presence 
of a source term a"(x), (Tj(x) or (Tij{x.) (or a time derivative thereof), the three-dimensional 
integral / c/^x extends only over the compact support of the material source. The finite part 
prescription is unnecessary for such terms. [Note that /|'"^ is identical to the 2PN expansion 
of the exact linearized gravity result given by Eq. (5.33) of [^.] The "F-terms" 



(named after the quantity defined in Eq. ( |3.22D below) denote all the contributions involving 
the three-dimensional integral of the product of (spatial derivatives of) two Newtonian-like 
potentials, e.g. 

FPb=o / d^^l^l^'xi.Ld^ldiUdjU] = ^ {fPb=o / d'^l^fxijLdiUdjU^ . (3.2) 

The time derivatives appearing in the F-terms have all been written in a manner such 
that they can be factorized in front as total time derivatives acting on (the finite part of) 
a three-dimensional integral. Finally, denotes the only term of ( |2.17| ) involving the 
three-dimensional integral of a term trilinear in source variables, namely 

4"^^ = -^FPb=o I d''^\^\^XLW,,di^U , (3.3) 



where we recall that Wij, defined by (|2.14c| ), is a bilinear functional of (T(x',t). 

We shall consider in turn the three contributions to II. The "Compact" and "F" contri- 
butions will be evaluated in the quadrupole case (£ = 2) while the "W^" contribution will be 
calculated for any £. The cases £ = and 1=1 play a special role as they do not correspond 
to radiative moments, but to conserved quantities. We check in Appendix ^ the agreement 
with known results for these low moments. 
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A. The compact terms and their expUcit form in the circular two-body case 



The general, A^-extended-body expression for the "compact" contributions to the 2PN 
mass moments reads 



4(2£+ 1) 



(£ + l)(2£ + 3)c2 
2(2£+ 1) 



XiL 



4f/ 



+ —^ss + 



1 /x2„, 1 X^ 



2(2£ + 3) Vc^^*'^^ 4(2^ + 5) c^"* 



^ , , 2U, 2Ui 



1 



2(2£ + 5) 6 



2 ^t^i 



+ (£+!)(/+ 2)(2£ + 5)c4^*^'^^'''*^; 



(3.4) 



From Eqs. ( |2.19| ) we obtain the corresponding point-mass form 



AC) 



A=l 

1 



2 



1 



(/^A(yA)'yi' 



2(2£ + 3)c2 

4(2^ + 1) d 
■(£ + l)(2£ + 3)c2rft 
2(2£+ 1) 



8(2£ + 3)(2^ + 5)c4dt4 
2t/^ 




1/r 



(£ + l)(2£ + 3)(2£ + 5)c4rft3 

2(2£ + l) rf2 
(£+!)(£ + 2) (2£ + 5)c4rft2 



{fiAV\v^AyA^) 



(3.5) 



in which we have introduced for brevity fiA = /iyi(l + v^/c^)- As it is written, the result 
( p.5| ) depends (at 2PN order) not only on the positions and velocities of the N 
compact bodies, but also on higher time derivatives thereof, up to d^jA/dt^. To reduce 
the functional dependence of to a dependence on positions and velocities only, we need 
to use the post-Newtonian-expanded equations of motion of the A^-body system. At this 
juncture, we restrict ourselves to the simplest case of a binary system evolving on a quasi- 
circular orbit. Consistently with the 2PN accuracy of the generation formalism considered 
here, we use the 2PN truncation of the 2.5PN equations of motion of Ref. [|lT|. We use 



a (harmonic) coordinate system in which the 2PN center of mass is at rest, at the origin. 
In such a coordinate system the 2PN truncation of the equations of motion admit exact 
circular periodic orbits. Using the 2PN-accurate center of mass theorem of Ref. |^ , we can 
express the individual center-of-mass-frame positions of the two bodies in circular orbits in 
terms of the relative position 



as follows: 



yi 
y2 



yi -y2 



[X2 + 3z/7'(^i-^2)]x, 
[-Xi + 3z/7'(^i-^2)]x. 



(3.6) 



(3.7a) 
(3.7b) 
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(Eqs. ( 
in Ref. 



and 



are obtained by setting Gi 
.) Here, we have denoted 

m = mi + m2 , Xi = 



for circular orbits, where Gi is given by Eq. (19) 



nil 
m 



7 = 



X1X2 = 
Gm 



77117712 



m2 
m 



= 1-Xi 



9 ' 
c r 



(3.8) 

(3.9) 
(3.10) 



with r 



|yi 



y2| denoting the constant (harmonic) coordinate radius of the relative 
orbit. One should note in Eq. ( p.7|) the absence of IPN corrections to the usual center-of- 
mass expressions. This is an accident due to the fact that we are restricting our attention 
to circular orbits. [In the non-circular case there are IPN corrections to X2 and —Xi which 
are proportional to — Gm/r, see e.g. ||2^.] Then the content of the 2PN equations of 



motion reduces to the knowledge of the 2PN-accurate orbital frequency a;2PN given by 

Gm 



2 

'^2PN 



(3 - z/)7 + 



41 
4 



7 



which is such that 



~dt 
dw 

Itt 



-^^pnX + 0(£^ 



(3.11) 

(3.12a) 
(3.12b) 



Let us note that Eqs. ( p.l2| ) imply, as usual, that t> = |v| = uj2PNr + 0{e^), so that Eq. ( p.ll| ) 
imphes 



7^ = ^ 



1- (3-z/)7 + (6 + ^z/ + z/2 



7 



(3.13) 



and we restrict 



(C) 

We are now in a position to compute explicitly the "compact" terms J)^ , 
ourselves to the quadrupole case i = 2. 

(c) 

Without entering into the details of the calculation of I^j , let us mention that there arise 
many symmetric functions of the two masses which can be straightforwardly expressed in 
terms of the total mass m and of the quantity u = X1X2 by using the well-known fact that 
a symmetric polynomial in Xi and X2 can be written in terms of the elementary symmetric 
combinations Xi + X2{= 1) and X1X2. Useful formulas for this reduction are: 



X^ + X^ 

x! + xl 
xt + xt 
xl + xl 



2v , 
3z/ , 

4z/ + 2z/2 
5z/ + 5z/^ 



(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 



The computation of the quadrupole 11^ for two bodies and circular orbits is long and tedious 
but quite straightforward. We obtain 
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I^P = STFij urn 



1512 



x^-' (5203 - 182751^ - 2785i^^) 



7 r 



t;^^(191 - 577Z/+ 109Z/2 



(3.15) 



where STFjj denotes the STF projection with respect to the indices ij. 



B. The quadrat ically nonlinear terms 

Separating out from Eq. (p.l7|) the terms involving the three-dimensional integral of a 
product of (spatial derivatives of) two Newtonian potentials, we define 

~d,U,djU, + 2d,U,d,U, - ^d^{U^)] 

-AidjUdiUj - -diUdjUj 



(£ + l)(2£ + 3)7rGc4 ' V ' ' 4 

AidiUdjU)] . (3.16) 



(£+l)(£ + 2)(2£ + 5)27rGc4 

These terms can be deduced from the results of Ref. [If. Denoting ri = |x— yi|, r2 = |x— y2|, 
= |yi — y2|, the latter reference has introduced the kernel 

A;(x; yi, y^) = \ ln[(ri + r^f - r^J (3.17) 

and proven that it satisfies (in the sense of distribution theory) the identity 

1 



A^A;-27r5i2, (3.18) 
rir2 

where 5i2 denotes a distribution supported on the segment joining yi to y2: 

^i2(x;yi,y2) = ri2 / rfa(5(x - y^,) , (3.19a) 

JO 

where 

yc, = (1 - a)yi + ay2 . (3.19b) 

Moreover, the kernel is such that its multipolar expansion when |x| oo contains, 
besides a logarithmic term, ln(2|x|), only terms of the type |x|~^~^P?7,i (with p G IN). The 
important point, as we are going to see, is that the multipolar expansion of /c(x) contains no 
terms of the type |x|~^~^n2,, i.e. no homogeneous solutions of the Laplace equation. [This 
is the feature defining the kernel k by contrast with the kernel g satisfying = ('"i'"2)~^ 
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everywhere and containing a homogeneous piece whose distributional source is precisely 



the — 27r5i2 term in Eq. (|3.18|) ; see [§.] The latter property implies that 

FPb=o / d^^\^fxLA^k{^) = FFb=o I c/'xA,.(|x|^xl)A;(x) 

= FPij=o 1^(5 + 2i+l) J rf^x|x|^+^-2niA;(x)| = . (3.20) 

Here, the first equality is obtained by integrating by parts (the surface term at infinity 
vanishing by analytic continuation from the case where Re{B) is large and negative), and 
the last follows from the fact that ^(x) is well behaved at the origin |x| = and contains no 
"homogeneous" terms \:s.\~^~^nL in its multipolar expansion at infinity. Indeed, going back 
to the definition recalled above of analytically continued integrals, as the sum of one integral 
over the ball < |x| < tq and one over its complement |x| > tq, we see that the factor B in 
front will give a zero result, except if the integral near the origin (0 < |x| < e) or the one 
near infinity (|x| > tq, with tq arbitrarily large) generates a pole at i? = in the complex B 
plane. As A;(x) exhibits no power-like blow up near the origin x = (even when the latter 
coincides with yi or y2) the integral near the origin is easily seen not to generate any pole at 
-6 = 0. Concerning the integral near infinity, written as dr r^~^^{Jg^ dQ hik), we see from 
the orthogonality of the his over the sphere 5*2 that only the presence of terms oc fiL/r^^^ in 
the multipolar expansion of /c(x) could generate a pole through ^B+e^-e-i _ 

By combining the identity (|3.18|) with the result (|3.20| ) we conclude that 

FFb=o I c/=^x|x|^XL^ = -27rr^(yi,y2) , (3.21) 
J rir2 

(see also Sect. IV in Ref. ||T0[) in which, following the notation of [Q, we have introduced 
F^(yi,y2) = J c/=^xxL5i2(x;yi,y2) = ri2^'rf«y<^> , (3.22) 



where Ua^^ denotes the STF projection of . . . Ua with defined by Eq. ( p.l9b|) . As the 
dependence of Ua^^ on a is polynomial it is easy to perform the integration over a in (|3.22|) 



to get m 



vLr \ yi - y2 {L-p p) oQ^ 
Y (yi,y2) = 2^yi 2/2 , (3.23) 



where 1/2 = 1/2^ ■ ■ ■ ifl 1 Ui'^ = Vi*^ ■ ■ - Ui ■ By taking derivatives of both sides of ( |3.21| ) with 
respect to y\ or ^"^^ integrating over yi and y2 after having weighted the integrand with 
some source functions ^^(yi), c'"/3(y2) [where cTq denote a, cTj or (jjj], one can obtain all the 



terms in (|3.16|) since these are all bilinear in some Newtonian-like potentials. For instance, 




the right-hand side of Eq. ( |3.2| ) can be written as 

c^Vi^(yi,t) / rfV2a(y2,t) {--^^ >^^^'^(yi,y2)| , (3.24) 

where one must be careful about the minus signs appearing in the spatial derivatives due to 
dr^^ jdx^ = —dri^/dy{, 9r^^/9x* = —dr2^/dy\^ and about keeping the total time deriva- 
tives factorized in front of the whole expression. It is convenient to introduce a special 
notation for the derivatives of with respect to yi and y2, say 
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d d 



' dy\ dyi 



Yh 



dyi dyi 



(3.25a) 
(3.25b) 
(3.25c) 



With this notation in hand, it is easy, from the resuh ( p.21| ), to obtain the following expres- 
sion for the "y-type" contribution to II defined in Eq. ( p.l6|) : 



r(^) 



■ sk 



4(2£+ 1) 



(^+l)(2£ + 3) 
2£+ 1 



4 



0"l 0-2 



Y 



aL 



2(£ + l)(£ + 2)(2£ + 5) 



abL 



(3.26) 



where cTi = cr(yi, t), (72 = o-(y2, t), etc. .. . 

Finally, the point-mass limit is obtained by inserting Eqs. (|2.19|) into Eq. ( p.26[) . Note 
that, because of the overall factor, it is enough to use the Newtonian approximation for 
the source terms [e.g. cr(x, t) = S^mA5(x — yA{t)) + 0{e^)]. To increase the readability of 
the result, it is convenient to introduce a shorthand for the contractions of the derivatives 
( p.25| ) with the velocities: 



■ VAVC 

— 7)* 1!^ ■ ■ 



where A, C 



1, 2, as well mixed notation like e.g. 



Y^ 



(3.27a) 
(3.27b) 
(3.27c) 



(3.27d) 



This leads, for any i, to the following expression for the F-terms 



r(^) 



^E^Amc{2n^,.,-i; 

A,c ^ 
1 + 2 Y^ - -(9^ (y^) 

r) VA^VC ' ^ VC^VA 9 *V J 



L 

VAVA 



4(2£+ 1) 



(£ + l)(2£ + 3) 
2i + l 



vc^ s /^^^ s 



sL 



2(£+l)(£ + 2)(2£ + 5) 



aYb 



abL 



(3.28) 



in which all functions Y^ are evaluated with in their first argument, and y^ in their 



second one, e.g. 



VA Vc 



L 

VA^A 



t>^f;^Y^^(y^, yc), and in which all the 
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self-terms A = C must be omitted. Finally, by using the expression ( p.23|) for Y^, a long 
calculation yields the following explicit expression for the y-type mass quadrupole for two 
bodies in a circular orbit, and considered in the center-of-mass frame: 



2mi"-f 
63 



STF 



7x^^ (55 - 155z/ - 53u^) + —v'^{-118 + 92z/ - lOz/^) 



(3.29) 



C. The cubically nonlinear term 

Let us now tackle the cubically nonlinear term 

= ^FPb=o / rf=^x|x|^x^W^,,9,,t/ , (3.30) 

with Wij = A^^{diUdjU). We shall show how to evaluate explicitly, for all values of i, 
in the case where the source is a binary system. [Note in passing that as ( |3.3CI| ) depends only 
on the instantaneous mass distribution of the source our result is valid whatever be the orbit 
(circular or not) of the binary.] To do this for all values of i we need what is essentially a 
generalization of the method of |p, i.e. a detailed study of some cubically nonlinear kernel. 
In the particular (and most urgently needed) case of the mass quadrupole, i = 2, one can 
evaluate ( |3.3CI| ) by other means, as is shown in Appendix ^ which succeeds in computing 
/f]^^ for N-body systems. This gives us an independent check of the results below. 



Let us define a function W{'x,t) [which is a trilinear, nonlocal functional of a{x.',t)] by 

AW = -^W,,d,,U , (3.31a) 

iy(x) = / -^W^.,(x')a.,f/(x') , (3.31b) 

where, for brevity, we suppress the dependence on the time variable which is the same 
on both sides. Note that the integral ( |3.31b| ) is a usual, convergent integral [at infinity 
Wij = 0(l/r), dijU = 0(l/r^)]. As usual, W is characterized as being the unique solution 
(in the sense of distribution theory) of ( pj.31a| ) which falls off at spatial infinity. 



Our method for computing I^^^ is similar to the one we used in Eqs. ( |3.17| )-( |3^20| ) above 
to compute Ij^^ from the results of |^ on the kernels k and g. Inserting ( |3.31a| ) into ( p.30| ) 
gives 

= -^^^B=o I rf3x|x|^x^A.l^(x) . (3.32) 

Integrating by parts the two spatial derivatives of A^, (the analytic continuation from a case 
where Re{B) is large and negative ensuring the vanishing of any surface term at infinity), 
using the formula A^(|x|'^x'^) = B{B + 2i + l)r^~^^~'^n^ , where r = |x|, n* = x'^/r, and 
writing explicitly the definition of the analytically continued integral (in polar coordinates: 
(i^x = r'^drdVt) we get 
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YVb=Ab{B + 21+1) 



roc r 

+ dr dVtr^+^fi^W 

Jrn J 



(3.33) 



Because of the regularity of Vr(x) near the origin x = the first integral on the right- 
hand side of ( |3.33|) will continuously depend on i? G (T near B = 0, and will therefore not 
contribute to because of the explicit B factor in front. We are therefore left with (since 
the second integral can have at most a simple pole) 



2/ -I- 1 ( r°° 

FFb=o<B I drr 



B+l 



dVL h^W{-K) 



(3.34) 



where Tq can be taken arbitrarily large, so that the expression ( p.34| ) depends only on the 
asymptotic expansion of iy(x) for |x| — oo, say 



iy(x)= y: 



p>l,f>0 



n 
rv 



(3.35) 



Although this is not a priori evident from its definition ( p.31| ), the explicit expression we shall 
derive below for iy(x) shows that the multipolar expansion ( |3.35| ) proceeds according to the 
inverse powers of r, without involving logarithms. [Actually, the presence of logarithms in 
( p.35| ) would not change the result below.] As was already mentioned above, the only terms 
in the multipolar expansion ( |3.35|) which can generate a simple pole in B so as to cancel 
the factor B in front are those which correspond to an homogeneous solution of the Laplace 
equation oc jr^'^^. Let us then define the "homogeneous" piece of W by restricting the 
double sum in ( p.35|) to the case where p = £ + 1, say (with the definition = W^l"*"^) 



W 



hom / 



£>0 



hom 



n 



(3.36) 



Note that this definition is independent of the choice of the origin x = around which one is 
expanding. This stems from the easily verified fact that under a constant shift x — >■ x + c the 



homogeneous and inhomogeneous pieces of W do not mix. Using 
-^-1 + 0(5°) and/dl^ninL/fi/ = [47r£!/(2£+ 1)!!]^^ (where n\\ 

r(W) 



ro ' 

n(n 



-B 



(n-2)(n-4)...(l or 

2)) we get the following link between I'^^ ' and the coefficients appearing in the multipolar 
expansion of iy^°'"(x). 



GI 



(W) 



(2£- 1)!! 



w 



hom 



(3.37) 



Inserting back (p.37 ) into ( 3.36|) shows that the multipolar expansion of ly'^"™ can be com- 
pactly written in a form which is familiar : 



W 



hom 



r{W) 



a'. 



(3.38) 
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The form ( 3.38|) is related to another, more convenient, expression for /|^'' [in terms of the 
"source" of ly"^"™] that we shall derive below. Let us now obtain an explicit expression for 
iy(x) in the point-mass limit and extract its homogeneous part. 

Inserting U = Gnii/ri + Gm2/r2 in the definition of Wij yields (with n\ = (x* — y{)/ri, 
etc.) 



AW,. = 



2n\ni 



i J 



rir2 



whose solution can be written as 



ij ] 



with 



11 



w., 



22 
12 



^5.1nri + ^^, 
Inr + 

i9j ~l~ j9i ■ 



(3.39) 



(3.40) 



(3.41a) 

(3.41b) 
(3.41c) 



In Eqs. ( ^.41aD , ( |3.41bD the derivatives = d/dx\ while in Eq. ( ^.41c| ) g is the quadratic 
kernel [HI 



5((x; yi, ya) = ln(ri + rs + r^) , 

^xg = — , 
rir2 

and we have introduced the same abbreviated notation as above: 

d d 



i9j 



dyl dyi 
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From ( |3.4CI| ) we get the cubically nonlinear effective source term 

Wi.dijU = G'imlA''' + mlm2A''^ + m.mlA'^' + mlA^^^ 

with 

-ri 8 - 
1 1 



11 11 

W^^d = -d-\nr,d-— = — 



A''' = w!fd., 



r2 



2igjdij h -dij Inridij h ^^A— 

Ti 8 r2 8r( r2 



(3.42a) 
(3.42b) 



(3.43) 

(3.44) 

(3.45) 
(3.46) 



the other terms in ( p.44| ) being obtained by exchanging the roles of yi and y2. 

Note that the two x-derivatives acting on 1/ri and l/r2 in ( p.46|) [which introduce non 
locally integrable singularities in x-space] are left uneffected and must be interpreted in the 
sense of distribution theory (in x-space). Rigorously speaking, one is not allowed to work 
within the framework of distribution theory (because one is dealing with the product of a 
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distribution dijr^ ^ by a function igj which is not smooth at x = yi). One should, e.g., use 



the well-defined analytic continuation procedure of |^ which has been shown to correctly 
describe the self-gravity effects of compact bodies. [The latter analytic continuation proce- 
dure yields in particular an unambiguous treatment of the self source term A^^^, Eq. (|3.45| ).] 
In practice, a technically easier way to deal with this subtlety is to work with the yi— and 
y2— derivatives of quantities which are less singular in x-space (see the first terms in the 
definitions of H and K below). When doing so, there arises only one term which is not well- 
defined in the sense of distribution theory, and this term, oc ni5(x— yi), clearly vanishes when 
treated more properly by analytic continuation. Computing A~^{WijdijU) is easy for A^^^ 
(using Arf'^ = 6rf^) and the last term in A^^"^ (using rf^Ar^^ = r{2^Ar^^ = A(rf2^r^^)). 
The other contributions to A^^"^ are much more intricate to deal with. We succeeded in 
evaluating explicitly A~^A^^'^ by combining the results of Refs. |2^j3^ which pointed out 
the usefulness of considering certain yi— and y2— derivatives of combinations of g, Inri, ri, 
r2 and ri2, with the fact (contained in a somewhat roundabout way in Ref. [^) that the 
inverse Laplacian of rf ^r^^ has a simple expression: 



A" 



rfr2 



1 r2 
2rf2 rf 



(3.47) 



[The latter being most simply obtained from its easily verified "inverse": A(r2/r^) = 
2r^2/(^i^2)-] Denoting Vi = d/dyi, Ai = d/dy{d/dy\, Vi-V2 = d/dy\d/dy\, and 
n\2 = iVi ~ yl)/^^^, we define the quantities 



if = Ai(VrV2) 

1 1 

+ o 2 + 

2 ri2rf 
3ni2-ni 



2 

ni2- ni 

2 ' 

run 



1 r2 



2 r?2r? 



+ 



K = 



+ 



ni2 „ 
^- Vi^ , 

— Inri lnri2 

r2 r2 

11 11 

2 r'i2r'i 2 r2 2 rf2r2 



2 rl2ri 
(VrV2) 
1 r2 



which verify 



= 2 j^fj dij — 



1 



A^K = 2diAnridij— . 



r2 



Finally, introducing the combination 



4 2rf2r2 

we find that W{x.), the unique solution of (|3.31a|) falling off at infinity, is given by 



(3.48) 
(3.49) 

(3.50) 
(3.51) 

(3.52) 
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w 



91 



— \ — m\m2Q + (two other terms obtained by exchanging 1^2) 
3ri 



(3.53) 



To project out from W{^) the part 1^*^°™ whose muhipolar expansion, when |x| oo, 
is purely homogeneous, it is convenient, besides using the defining criterion that it contains 
only terms of the type n^/|x|^+^, to notice that 1/1/'^°™ must also be a non smooth function 
of yi and y2 (considered jointly). This can be shown either from the general structure of our 
generation formalism, or, more simply, by remarking, on dimensional and tensorial grounds, 
that I'l^^ must be of the form: G'^m\m\~^ / {c^r\2) times a tensor product of (. vectors yi or 
y2. When using these two criteria in conjunction, one finds that many terms in H and K 
project out to zero. For instance, by explicating 



Ai 



-9 



1 

ri 



1 



9 



1 

ri 



1 

ri2 



r2 



(3.54) 



with the decomposition g = k + \h (see Ref 

^(x; yi,y2) = In 



^2/ ri ri2 ri2ri 
^.u y^^^ ^) where k is smooth in (yi,y2) and where 
ri+r2 + ri2 



ri + r2 
47r(5i2 , 



ri2 



(3.55a) 
(3.55b) 



has a purely homogeneous multipolar expansion, one finds that the first term in H 
(Eq. ( |3.48|) ) projects out to zero. Finally one gets 



W 



horn 



15 / r2 



horn 



ni2 ni 



8^2 



12 



5 1 



8 rf2r2 



^"12 



ri2 rf 

+ (1^2) 



(3.56) 



where (r2/r^)^°™ denotes the homogeneous projection of r2/r^ and where h was defined in 
Eqs. ( p.55| ). Let us now evaluate u = (r2/r^)'^°™. By expanding r2 = |x — y2| = |ri + ri2| 
[with ri = x — yi, ri2 = yi — y2] in powers of ri2, we get the expansion at infinity of r2/r\: 



!l - V — 

2 ~ £1^" 

' 1 i>o ^■ 



(3.57) 



[where, as usual, rfg = Ti2 • • • "'^12]- The projection of dLTi/rf onto the homogeneous solutions 
oc nf /r{'^^ (using x = yi as origin for the expansion at infinity) is simply obtained by taking 
the STF projection of the multi-index L. Hence (with (—1)!! = 1, (—3)!! = —1) 



r2 ' 



horn 



2^ /I ^12 ^2 



r 



-(2^-3)!!r 



"1 



(3.58) 



Using formula (A25) from [^ this can be rewritten in terms of Legendre polynomials of 
nini2: 
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« = — E y^^^K- ^12) i—T' . (3.59) 



It is clear from Eqs. ( 3.58|) - (|3.59|) that u(x, yi,y2) is a solution of A^m = which has an 



axial symmetry around the straight line joining yi to y2. We can represent u in closed form 
by introducing a distribution of "charges" along the segment yi — y2, say 

/■I wia) 

M(x;yi,y2)=/ da- -, y^ = (1 - a)yi + ay2 • (3.60) 

Jo |x-ya| 

By identifying ( p.60| ) and ( p.59| ) on the axis of symmetry, we see that the weight w{a) with 
which the "charges" are distributed on the segment yi — y2 must satisfy 

1 X 

, daw{a)a^ = —— — -. (3.61) 
2t — 1 

We find that w{a) must be defined within the framework of distribution theory (rather than 
that of ordinary, locally integrable functions) as 



w{a) = Pf 



(3.62) 



where the symbol Pf denotes Hadamard's partie finie for an integral over the interval 
< a < 1 (see Ref. ||3^). Actually, this could also be written in terms of the analytic- 



continuation Finite Part operator used in our general formalism, but Hadamard's P f oper- 
ator is a simpler object when, as is the case in Eq. ( p.62| ), the integrals to be defined involve 
non-integer powers. We want also to emphasize by the change of notation that the appear- 
ance of the distribution w{a) is quite disconnected from the formalism behind our starting 
formula ( |2.17| ). We note in passing that we could dispense of using Hadamard's partie finie 
by using w{a) = {d/da){a~^^'^), where the derivative is taken in the sense of distributions, 
and by integrating ( |3.6U| ) by parts. However, the form ( |3.62| ) is technically more convenient 
for our purpose. 

Summarizing the results so far, we have succeeded in resumming the infinite multipolar 
series defining 1^^^°™, Eq. (|3.36| ), to obtain it as a finite sum of derivatives of rf^ and r^^ [see 
Eq. (|3.56|) ], plus two more complicated homogeneous solutions: a yi-gradient of the function 
h, Eqs. ( |3.55| ), and the function u, which is a distributional superposition of elementary 
solutions |x — yal"^ on the segment yi — y2. Note that while Eq. ( |3.36|) defined only the 
asymptotic expansion at infinity of ly*^"™, our closed-form result for W^°'^{x.) defines it for 
all values of x. At this stage, it is convenient to bypass the link ( p.37|) - (p.38|) between the 
object we seek, and the coefficients of the multipolar expansion of iy'^°™(x) by using 

our detailed results on ly''"™ for defining the distributional "source" (in x-space) of 1^^^°™ 
by 

^(x) = -i-A,W^'^°-(x) . (3.63) 

Thanks to our resummation of TV*^"™ as a sum of elementary "homogeneous" (in the sense 
of functions, but not of distributions) solutions, the definition ( ^.631 ) makes sense with S'(x) 
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being a distribution (in x-space) whose support is localized on the segment joining yi to y2. 
More precisely, we have 



S 



m\m2 



15 1 



-Vr V2(ln ri25(x - y2)) 



i« 2-Pf 
16rf2 

"12 



da a-'^l^5{^ - Ya) - 4—- Vi(5(x - yi 
Jo ri2 



— ^o(x 



y2j 



, ni2 



Vi ^ri2 j da 5(x 



2) 



(3.64) 



The introduction of the source S of 1^'^°™ 

t{W) 

deduce as usual 



simplifies our evaluation of the multipole moments 
ly ' . Indeed, from the definition (|3.36|) and the fact that 5'(x) has a compact support we 
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By comparing with ( p.38| ), we find 
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(3.66) 



It is interesting to remark that the work of the present subsection, together with the one 
of the previous subsection, is analogous to the method used in Ref. p[. Basically, we have 
decomposed —4WijdijU/c'^, which was an effective source term for the inner metric (i.e. a 
cubically nonlinear analog of the quadratically nonlinear terms obtained from the elementary 
object rf ^r^^) into a piece A^W^"^ = —AtiS which has a compact support (analog to A^h = 
—4:7t6i2, in 1^, defining the compact source rj^'^) and a complementary piece A^W^^^"'^ , with 
y^mhom ^ _ yf/hom i^gj^g i\^q aualog of the kernel k, which does not contribute to ( p.lTI) . 
In other words, the combination of the results of Ref. [|TU| with the results given here, end up 
by yielding a formula for given explicitly by an integral over a compact support (after the 
two replacements rf^r2"^ — > — 211612 and WijdijU / c"^ — > nS in Eq. ( [2.17] )). From Eq. ( ^.66[ ) 
one can read off directly from ( p.64|) an explicit expression for Note that one can get 

rid of the logarithm appearing in the first term of this expression by explicating the Vi 
derivative: 



-(VrV2)(lnri2j/2'') = 
This leads to our (essentially) final result 
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where, as usual, = y^^^ ■ ■ -y^^, and where, because of the dissymmetry betweeen yi and 
y2 in the last integral, it is important to recall the definition used here: Yq, = (1 — Q;)yi + ay2 
(which differs from the notation used in ^ by the replacement a — 1 — a). 

To get a completely explicit expression for Ii^^ in terms of yi and y2 one needs first to 
compute the integrals appearing in (|3.68|) : the first one has been given in Eq. ( p. 231) above 



(actually, one uses here the equivalent expression given by Eq. (5.24) in [^), the second 
one is simply obtained by expanding the tensorial power y^ = (yi — ayi2)'^^ in powers of 
a and by using the elementary integrals ( p.61| ). Then the derivatives with respect to yi 
and y2 must be effected. One can drastically cut down the whole calculation by using the 
fact that I^"^ must (as is clear from ( |3.66| )) transform under global spatial shifts of the 
coordinate system, x ^ x + e, y^ — * y^ + e according to 6sl^^ = ^£<ii,lL^l^ + 0{e'^) (see 
also Appendix B of [§). Defining for brevity some ^-dependent coefficients, 

Q = 25° + 45l + (-)^^^±^ (3.69a) 

(where 5^ = 1 if £ = and otherwise, and similarly for 6j), we obtain finally 



p^o\pJ \pj pK^-pV- 



(3.69b) 



where the multi- indices are such that there are i — p indices on yi and p on yi2 = yi — y2- 
Writing out explicitly the low orders in i we find 

'r = + (1 " 2) , (3.70) 



with 



'12 



Q = 1 , (3.71a) 

Q^ = yl + yi2, (3.71b) 

Q,,. = yf'> + 2yfyi> + 2yf^'> , (3.71c) 

Q^,k = yf'"^ + ^yf'y'l^ + Qyfyi'^ - 22/it^'> , (3.7id) 



Q^JM = yf'"'^ + "^yf^'yil + i^yf'y'A^ - Syfyi'^'"- + ^-^yT'"- ■ (3.71e) 



7 

One should note that in all the formulas for W and I^^'^ given above, starting with Eq. ( p. 531 ), 
one must complement the explicitly written results by adding similar terms, proportional to 
mim\, obtained by exchanging the labels 1 and 2 (remember that y^i = y^, ~ = "~Z/i2)- 
For instance, the complete expression of the quadrupole reads 

) = - ^Yr' STF., [m,y^ + m,y^ + 2{m,y\ - m^y^^y{^ + 2{m, + m^)y%] . (3.72) 

An independent, direct derivation of the quadrupole ( |3.72| ) (without using the decomposition 
of W into homogeneous and inhomogeneous parts) is given in Appendix 0. In the center of 
mass frame, we finally get for the cubically nonlinear contribution to the quadrupole 
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/^^) = -z/m72(2 + hu)x'^ . (3.73) 

Summing up the explicit results obtained for the three pieces of the 2PN-accurate quadrupole 
(considered for circular orbits, in the center-of-mass frame) we get 

2 ij 

- 1:^(461 + 18395z/ + 241z/2) 
1512 



2 ij \ 

+ ^^(1607- 1681z/ + 229z/') . (3.74) 
378c^ 



IV. THE 2PN-ACCURATE WAVEFORM AND ENERGY LOSS 



A. The waveform including its tail contribution 



The 2PN-accurate waveform is given by Eq. (|2.3|) in terms of the "radiative" multipole 
moments Ul and Vl which are in turn linked to the source moments II and Jl by Eqs. ( p.6|) 
and ( p.7|) . The latter equations involve some tail integrals and therefore yield a natural 
decomposition of the waveform into two pieces, one which depends on the state of the binary 
at the retarded instant Tr = T — R/c only (we qualify this piece as "instantaneous"), and 
one which is a priori sensitive to the binary's dynamics at all previous instants Tb, — t<Tr 
(we refer to this piece as the "tail" contribution). More precisely, we decompose 
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where the "instantaneous" contribution is defined by 
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(4.2) 



and where the "tail" contribution reads 
2G 2Gm 



(^fcm)tail 
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3c \2b2 
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+ -\n( — )Njl^>iTn 



+Yc 1^ (i) ^'^K^N,J^?aiTR - r)} . (4.3a) 



We have used for simplicity the notation 
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Note that the decomposition ( [4.1| )-( ^y3D into instantaneous and tail contributions is conve- 
nient but by no means unique. In particular the logarithms of the constants bi, 62, &3 in 
Eq. ( [4.3a| ) could be as well transferred to the instantaneous side ( |4.2| ). 

We first study the instantaneous contribution ([4.2| ), which is straightforwardly computed 
from the 2PN- accurate mass quadrupole moment previously derived (Eq. ( pj.74D ), and from 
the knowledge of the other moments necessitating the IPN accuracy at most. The IPN- 
accurate source moments have been derived in ^ (see Eq. ( |2.9| ) above) and Note that 
the result (5.18) of ^ is equivalent to inserting Eq. ( |3.21| ) into the result ( |2.13| ) above. We 
list below the relevant mass-type and current-type moments which have to be inserted into 
the waveform ([4.2|) : 

■y •• 11 

^ (1 + 39z/) x'' + 

42 ^ ^ 21 c2 



= urn STF,, <j x'^ + 39z/) x'^ + — -v'' (1 - 3z/) 



2 ij 

-^-^ (461 + 18395Z/ + 2411/2 

1512 V 

+ (l607- 1681z/ + 229z/2) , (4.4a) 



= vm (X2 - Xi) STF,,, Y''' - 7^^^*'' + ^ - 2^^) j . (4-4b) 
= STF,,M [x'''\l - 3u) + ^x'^"' (3 - 125z/ + 345//^) 

+ ^ ^ v'^x^'il-^u + ^u^)] , (4.4c) 

I^Jklm = vm{l - 2z/)(X2 - Xi) x'^^'^ , (4.4d) 

I^Jktmn = ^^^(1 - bu + bu^) X*^'^'™" , (4.4e) 

J,, = z/m(X2 - Xi) STF,, e,ab [x'^'v" + ^ x^^'v^Gl - 8z/)| , (4.4f) 
Jijk = vm STFi.fc j x"*^t;^(l - 3z/) + ^-^(181 - 545z/ + 65//^) 
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+ 77 ^ x'^v'''^ (1 - 5z/ + 5z/2) !> , (4.4g) 
45 & J 

Ji.-fc; = vm{\ - 2z/)(X2 - Xi) STF.^H^Zafex'^^^^^i;'' , (4.4h) 

JiiWm = z/m(l - 5z/ + 5z/2) STFijM^ e^^t a;"^-''^'t;'' . (4.4i) 

Note that Ref. P] gave two different forms of the current-type moments: a "central" form, 
its Eq. (5.18), equivalent to the form (2.13) above (taken from ||10[)) and a "potential" form. 
We found the first form simpler to evaluate when one is interested in computing the moments 
in the center-of-mass frame. [By using it, we detected an error in ||33| which computed the 
potential form of the octupole Jij^. the coefficient of the twelfth term in Eq. (9) of Ref. |33 



should read —5/2 instead of +3/2 (this twelfth term was the only one not checked by the 
transformation law of multipole moments under a constant shift of the spatial origin of 
the coordinates).] The relevant time derivatives of the moments (|4.4|) , as well as all their 
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contractions with the unit direction N which are needed in the computation of (^fcm)inst, are 
relegated to Appendix p. One gets for the 2PN-accurate instantaneous part of the waveform 
(in the case of two bodies moving on a circular orbit): 
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where the various post-Newtonian pieces are given by 
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(4.5b) 
(4.5c) 



(4.5d) 



(4.5e) 



20(nAr)-^^^(l-5z/ + 5z/^ 



(4.5f) 



Here n is a shorthand for x/r. Usual euclidean scalar products are denoted e.g. by [Nn) = 
N- n. We recall that the parameter 7 = Gm/{c^r) is related to the relative velocity v 
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of the bodies by Eq. (|3.13| ). When comparing the waveform (^4.5|) with the less accurate 
1.5PN waveform obtained in Ref. we find several discrepancies between the coefficients 
entering the 1.5PN piece ^fj'^^ (Eq. ( [4.5eD ). We have checked directly, by differentiating and 
squaring the complete 2PN waveform ( [4.5|) for circular orbits, that the correct instantaneous 
part of the energy loss is recovered (i.e. Eq. ( [4.12|) below). Therefore we suspect that the 
contentious coefficients in Ref. are erroneous. 

The tail contribution (4.3) is more difficult to evaluate than the instantaneous one be- 
cause it involves an integral extending over the whole past evolution (or past "history") 
of the binary. In another sense, however, it is easier to evaluate because it necessitates at 
the 2PN-level the knowledge of the quadrupole and octupole moments of the binary with 
Newtonian accuracy only. This is clear owing to the explicit powers of 1/c appearing in 
Eq. ( |4.3a| ). (Only at the 2.5PN level shall we need to include a post-Newtonian correction 



into the tail contribution.) In order to compute the tail integrals, we shall follow a procedure 
which is a priori dangerously formal (although most natural), but has been fully justified 



in Ref. where it was proved to yield the correct numerical value of the integrals, pro- 
vided that a weak assumption concerning the behavior of the gravitational field in the past 
(— r — oo) is satisfied. This assumption is essentially that the £th time derivative of a mo- 
ment of order d. tends to a constant when — r — >■ — oo. It serves to preclude, for instance, the 
emission of a strong burst of radiation in the remote past which would make a non-physical 
contribution to the tail integrals (see also Ref. for a discussion). This assumption is 
satisfied in the case where the binary is formed by capture of two bodies moving on an initial 
quasi-hyperbolic orbit with small enough energy. The method consists (i) in substituting 
into the tail integrals in Eq. (|4.3a|) the components of the moments calculated for a fixed 
(non decaying) circular orbit whose constant orbital frequency is equal to the current value 
of the frequency at time Tr, i.e. oj = W2pn = i^2Pn(^_r)- Then it consists (ii) in evaluating 



each resulting integral by means of the formula 
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dr In — e'^^ = — 

\2bJ n 



^sign{n) + i{\n{2\n\b) + C) 



(4.6) 



where C = 0.577- •■ is Euler's constant and Q is the frequency of the radiation (i.e. a 
real number Q = ±nuj; sign(i7) = ±1 and \Q\ = nuo). Note that this formula is to be 



applied even though the left-hand-side of ( |4.6| ) is an absolutely convergent integral only 
when Vt possesses a strictly positive imaginary part, Im{Q) > (compare the appendix A 
and Sect. 3 in Ref. [^). As proved in the appendix B of Ref. [^, the method (i)-(n) is 



valid for an orbit which is actually decaying, and thus for which ci;(Tr — r) tends formally to 
zero when — r —>■ —00. The numerical errors made in following (i)-(ii) have been shown to 
be of order 0[^ln,^], where ^ is the usual adiabatic small parameter describing the decay of 
the orbit [C, ~ uj/uj'^), taken at the current instant Tr, i.e. C, = ^{Tr). The proof consists in 
showing that an adequately defined "remote past" contribution to the tail integrals is itself 
of order 0[^(Tr)], so that only the "recent" values of the frequency, near the current value 
uj = uj{Tfi), contribute to the tail integrals. Although this has not rigorously been shown, 
the proof could in principle be extended to orbits which had a non negligible eccentricity in 
the past, and had in fact whatever behavior in the past which is consistent with the weak 
assumption made above. Note that the decay of the orbit is driven by radiation reaction 
effects which are of order 0{e^) in the post-Newtonian parameter e, and so the adiabatic 
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parameter ^ is itself of order 0{e^). One can therefore safely neglect in the 2PN waveform 
all the errors brought about by the above procedure (i)-(ii). 

We now calculate the two independent polarizations associated with the tail contribution 
( [4.3|) with respect to two unit directions P and Q perpendicular to N and such that N, P, 
Q forms a right-handed orthonormal triad. We adopt the usual convention that P and Q 
lie along the major and minor axes of the projection of the circular orbit on the plane of 
the sky, respectively, and we denote by c = cosi and s = sinz the cosine and sine of the 
angle between the line of sight N and the normal to the orbital plane (c = cosi is not 
to be confused with the speed of light which we denote exceptionally by Cq in Eqs. ( [4.7D - 
( [4.10| )). Denoting by (p the instantaneous orbital phase of the binary (defined as an angle, 
oriented in the sense of the motion, such that = f (mod 2it) when the relative direction 
of the two bodies is n = P), and using the relevant time-derivatives and contractions of 
Newtonian moments as readily calculated from Appendix |C|, we follow the items (i)-(ii) 
above and bring the two "plus" and "cross" polarizations (/i+)taii = \{.PiPj — QiQj) {hfj^)ta.ii 
and (/ix)taii = liPiQj + QiPj){hJ^)tia\ into the form 
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and 



ih.U = I dr[-8cln (-) sin[20(T, - r)] 



- - In cos[0(T/j - r)] cos[0(Tr - r)] 
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(4.8) 



Still following (i)-(ii) we now compute ( [4.7| )-( ^^ by replacing by a linear function of time, 
0(T) = uT + 00 where u = co'(Tr) is the current value of the orbital frequency, and by 
applying the formula (|4.(j|) to each integrals. As a result we get 
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The orbital frequency oo and phase in ([4.9|)-( ^?T0| ) denote the current values oo = uj{Tjij 
and = 4>(Tji). Both a; and will be computed in the next subsection (see Eqs. ([4.28|) - 
( [4.30| )). [Evidently, though during the computation of Eqs. (|4^)-( |^rT0| ), has been replaced, 
following Ref. p5|, by a linear function of time, its actual time variation is nonlinear.] 



B. The energy loss and the associated binary's orbital phasing 



The 2PN-accurate energy loss given by Eq. ( ^.5| ) is split, similarly to the waveform, into 
an "instantaneous" contribution and a "tail" one. Let us deal first with the instantaneous 
contribution, which is defined by 
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All the needed moments to compute this contribution have been given in Eqs. ([4.4|) . Their 
relevant time-derivatives and "squares" can be found in Appendix C. By a quite straight- 
forward computation we get 
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As for the tail contribution, it reads (when retaining only the terms which contribute to the 
2PN order) 
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Note that at the 2PN order there appears only the tail integral associated with the mass 
quadrupole moment Jjj, which can be replaced by its usual Newtonian expression. We 
now evaluate Eq. ( |4.13| ) by the same method (i)-(ii) as used previously for the tails in the 
waveform. Namely we replace in ( |4.13|) the third and fifth time-derivatives of the quadrupole 
moment by Newtonian quantities valid for circular orbits, and we perform explicitly the 
contractions between these moments, being careful that one of them is taken at the current 
instant Tr while the other is taken at the former instant Tr — r. This readily brings ([4.13 
into the form 
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where for this computation can be assumed to be a linear phase ||3^ . Then the use of the 
integration formula (|4.6|) immediately yields 
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There is no dependence on the constant b parametrizing the freedom in constructing the 
radiative coordinate system. This is not surprising because the constant b (or rather bi) 
enters a term which is a total derivative (see Eq. ( [4.13D ) and thus, as we aheady noticed 



T4 , which vanishes identically for circular orbits. 

A central result of this paper, namely the complete 2PN-accurate gravitational energy 
loss rate from a compact binary moving on a circular orbit, is now obtained by adding 



Eqs. (|4l^ ) and (glSl) : 



dT, 



R 



32c^ 
'~5G 



1-7 



/2927 



I \- -i^] + 47r7^/^ + 7- I 

V 336 4 / ' ' V 9072 



/293383 380 



+ 



9 



(4.16) 



The 7-parameter is 7 = Gm/{rc^) where r is the harmonic radial coordinate (see Eqs. (|3.10|) - 
( ^.13| )). It is to be noticed that in the form ( |4.16| ), i.e. when the post-Newtonian expansion 
is parametrized using 7, there is no term proportional to u"^ in the relative 2PN contribution. 
This fact is somewhat surprising because all separate pieces making up the energy loss, i.e. 
all the different "squares" of moments listed in Appendix C, do contain terms proportional to 
z/^. However the final coefficient of z/^ in the 2PN correction term turns out to be zero. The 
expression of the energy loss ( [4.16| ) completes several previous investigations having obtained 
either the lower-order FN corrections or the limiting case where the mass of one body is 
negligible as compared to the other one (limiting case u ^ 0). The OPN leading term in 



Eq. (|4.16|) was known from Ref. |3^, and the IPN correction was added in Refs. p7| , p8| . The 
1.5PN tail correction [Air term) was first computed in the limit z/ — > |[T^ and then shown to 
be also valid for arbitrary mass ratios in Refs. |]35| , |39| based on Ref. 0. [The generalization 



of the energy loss expression for non-circular (eccentric) orbits has been obtained in Ref. |5C 
for the OPN leading term and in Ref. for the IPN corrections. Note that for non-circular 
orbits the tail term 4717^/^ has simply to be multiplied by a function Lp(e) of the eccentricity 
of the orbit, but unfortunately this function probably does not admit a closed analytical 
form (see Ref. |^5|] for the numerical graph of this function).] Finally the 2PN correction 
term in Eq. ( 4.16 ) was known, up to now, only in the limiting case z/ — 0, where it was 
computed first numerically [|l5l and then analytically |jl6|. For comparison with the latter 



references, and for later convenience, let us re-express the energy loss ( |4.16| ) in terms of a 
new post-Newtonian parameter defined by x = (Gma;2PN/c^)^^^- The 7-parameter is related 
to the x-parameter by the inverse of Eq. ( |3.11| ), namely 
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(which also does not involve z/^-terms). Inserting ( [4.17| ) into ( [4. 16]) (and keeping consistently 
all terms up to the 2PN order) readily yields 
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This expression (which we notice involves now a z/^-term) relates the two coordinate- 
independent quantities dEs/dTfi and u = ui2pni and is therefore the same in all coordinate 
systems. It may thus be compared directly with the expression obtained in Refs. [0,|1^ in 
the limit z/ — > 0. We find that the coefficient —44711/9072 agrees with the latter references 



|2| 



Let us now denote the 2PN-accurate Bondi energy loss rate ( [4.18| ), or total gravitational 
luminosity, by Cb = —dEB/dTn, and let us equate, by a standard argument, the rate of 
decrease of the dynamical energy £ of the binary system to the opposite of Cb, i-e. 



dS 



dTi 



-C 



B 



(4.19) 



R 



We shall admit here the validity of the balance equation ( f4.19| ) to the 2PN order, although 
it has been validated only to the leading OPN order and in the case of the binary pulsar 
T^. [Note that the 2PN order is the last order at which the binary admits a conserved 
energy S (and a conserved angular momentum J').] The balance equation (|4.19|) drives the 
variations with time of the instantaneous orbital frequency u and phase of the (quasi- 
circular) decaying orbit. These uj and are the ones which enter the expression of the 
waveform (see e.g. Eqs. (|4.9| )-( ^?T0| )). To compute u and one must evidently specify first 
what is the left-hand-side of Eq. ( 4.19 ), i.e. one must know the 2PN-accurate expression 
of the center-of-mass energy S for a fixed (non-decaying) circular orbit as a function of the 
parameter x. 

Let us recall from Ref. (extending Ref. [^]) that the 2PN motion of a binary system 
moving on an eccentric orbit admits the following representation: 



n{t — to) = u 

r = ar{l — Cr cosm) , 



et smu + — sinv + —{v — u) 



F G 
K Iv ^ — 7 sin 2f H — - sin 3f 



(4.20a) 
(4.20b) 
(4.20c) 



where 



V = 2 arctan 




(4.20d) 



The time t, separation r between the two bodies and polar angle (or phase) in these equa- 
tions correspond to harmonic coordinates. The motion is parametrized by some "eccentric" 
and "real" anomalies u and f , and ten constants enter Eqs. ( 4.20p besides the initial to and 
0o: n = 271 / P where P is the time of return to the periastron (or period); K = 0/27r where 
G is the angle of return to the periastron {K — 1 is the relative periastron advance per rota- 
tion); the semi-major axis a^; three types of eccentricities et, and e^; and four constants /, 
g, F and G entering purely 2PN terms. Most important for our purpose are the expressions 
of the constants n and K obtained in Ref. (see Eqs. ( p.llD -( prT^ there) in terms of 
the constant center-of-mass energy S and angular momentum J'. Denoting S = umS and 
JT" = Gvrn^h we have 
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(4.21) 
(4.22) 



The other constants in Eqs. ( [4.20[ ) have been computed in Ref. |^^; however we shall here 
only need, in addition to Eqs. ( 4.21 ) and ( 4.22|) , the fact that in the circular ox\y\i case, which 
is defined by = 0, the two other eccentricities et and vanish to 2PN order, hence the 
two anomalies u and v agree to this order, and that the three terms involving the constants 



/, F and G also vanish (see [^). The 2PN-accurate circular motion of the binary is thus 
simply described by the equations 



r = , 

i)Q = nK{t - to) , 



(4.23a) 
(4.23b) 



which show that the orbital frequency u in the circular orbit case is equal to the product 
nK. Now uj = nK can straightforwardly be obtained from Eqs. ( 4.21 )-( 1.22 ) and the fact 



that £ and h are related to each other by the IPN-accurate relation 1 + 28^^ = |(9 + z/)£^/c^ 
(consequence of = 0). We find 
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1, ,S 3 / 11 9\ ^ 



(4.24) 



The inversion of Eq. (142^) then yields our desired relation between the energy S = vm£ 
and the parameter x = {Gmuj / c^Y^^ for circular orbits: 
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With Eqs. ( |4.18| ) and (|4.25|) in hand, it is now a simple matter to transform the energy 
balance equation ( [4. 19]) into the ordinary differential equation 
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where we have introduced for convenience the adimensional time variable 
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(4.26) 



(4.27) 



Solving ( |4.26| ) leads to the variation in time 6 of the parameter x and hence of the in- 
stantaneous frequency u of the quasi-circular orbit. Similarly, integrating d(f) = oodT^ = 
{b/v)x^f'^d6 leads to the variation of the instantaneous orbital phase cj). We finally obtain 



Gmuo 



-3/8 



1 + 



; 10 

37r 
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(4.29) 
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where 
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(4.30) 



and where 0c and 6c denote the values of the phase and adimensional time ( |4.27| ) at the 
instant of the coalescence. The frequency and phase ( [4.28D -( P:.29D have to be inserted into the 
expression of the waveform. Note that the coefficient of the 2PN contribution in Eq. (f4.29|) 
(proportional to q), increases by 52% between the test-mass limit (z/ = 0) and the equal- mass 
case (z/ = 1/4). This shows that finite mass effects (which cannot be obtained in perturbation 
calculations of black hole spacetimes) play a very significant role in the definition of 2PN- 
accurate theoretical waveforms. This proves the importance of post-Newtonian generation 
formalisms for constructing accurate templates to be used in matched filtering of the data 
from future gravitational-wave detectors. 



ACKNOWLEDGMENTS 

One of us (BRI) would like to acknowledge the hospitality of IHES during the initial 
phase of this collaboration. 

APPENDIX A: THE MASS MONOPOLE AND DIPOLE MOMENTS 



Since the expression ( |2.17| ) of the 2PN-accurate mass multipole moment is valid for all 
orders of multipolarity i, it is important to verify that for the lowest orders i = and i = 1 
it reduces to the expressions of the conserved mass monopole and mass dipole moments. 
This verification has already been done in Appendix B of [|l^] in the mass monopole case 
i = 0. Here we shall check that for two bodies moving on a circular orbit the mass monopole 
reduces to the expression computed in p6[, and the mass dipole is zero, as it must be since 
we are using a mass-centred frame. 

The mass monopole I and dipole /j, which are obtained by setting i = and i = 1 
in Eq. (|2.17|) , are split into three "Compact", "Y", and "W" contributions according to 
Eq. ( p.l| ) and evaluated separately using the method developed in Sect. |T1|. The "W" 
contributions have in fact already been computed in Eqs. (|3.70| ) and ( p.71a| )- (|3.71b|) . We 
quote only the results. For £ = we find 
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1 1 , 
1 z/7 z/(l 

2 ' 8 ^ 
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Adding together these contributions leads to 
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1 1/7 H — z/(7 — z/)7^ 

2 ' 8 ^ ^ . 
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which is in agreement with Eq. (20) in Ref. when specialized to circular orbits. For 
i = 1 we get 
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i 


= (X2 


r(y) 

i 


= (X2 


{W) 

i 


= (X2 



2 2Q \ • 

- + (A3a) 

(X2 -Xi)mz/72xS (A3c) 

in which we have used (for the "C" term only) the 2PN-accurate mass-centred frame equation 
( ^■71 ). The three contributions ( |A3D sum up to zero, 



/, = , (A4) 



as was to be checked. 



APPENDIX B: ALTERNATIVE DERIVATION OF THE VF-TERM 

We have computed in Sect. pTO of the main text the cubically nonhnear 'W term 



(defined by Eq. (|3.3CI| )) for all values of the order of multipolarity £, but in the special case 
where the source is a binary system. In this appendix we present an alternative derivation 
of this l^-term which is valid for a general fluid system (and hence for a system made of N 
compact bodies), but is limited to the quadrupole case £ = 2. When N = 2 and i = 2, i.e. 
when one evaluates the M^-term for a binary system and in the quadrupole case — this is 
what interests us in this paper — , we find that both derivations agree on the result given 
by Eq. 

Specializing the definition ( |3.3U| ) to £ = 2, we thus want to compute 



= ^7TFPij=o / d^^l^l^'x^.Wk^dk^U , (Bl) 



Gc^ 



where we recall that the potential Wkm is defined by 

1 r d^ 



W.mM = -— I -^MUdmU]{^,t) . (B2) 

47r J X — X 



We first perform on ( pi|) two integrations by parts in order to shift the spatial derivatives 
acting on U to the left side of the integrand. As can be proved thanks to Eq. (^.2|) of 



Ref. ||T^, all the terms coming from the differentiation of the analytic continuation factor 
|x|^ and having explicitly i? as a factor are zero. [Indeed, recall that Eq. ( |4.2| ) of Ref. |10 



permits one to freely integrate by parts all terms in the source moment ( p.l7| ) as if the 
analytic continuation factor and the "finite part" were absent.] The surface terms are also 
zero by analytic continuation. Hence we can split the Pi^-term into two pieces. 



1^ = 1^ + 1^^ (B3) 



pven respectively by 

2 



= ^ STF,, FPb=o / d^^\^fW,,U , (B4) 
4^^'^ = ^FPiJ=o / d^l^\^f[A^uWk<^x,> + dk^Wk^x,,]U (B5) 
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(where the angular brackets < > denote the tracefree projection). In ( p5|) the divergences 
of the potential Wij are given by 



dkW,k = -^^iU^) + - y j;^r^^H,U - f/5,a](x', t) , (B6a) 
9fc^l^fc„ = ^A(f/2) + -y ^— -^[aA[/-[/Aa](x',t) . (B6b) 

The most interesting contribution to compute is the first one, i.e. 1^^^^ given by (|B4D . 
This contribution has no explicit dependence on x in the integrand; this fact, which is 
special to the quadrupolar case, allows the alternative derivation followed in this appendix. 
To compute /j-]^^'* we employ the kernel g{'x; yi, y2) already introduced in Eqs. ( ^.42| ), 

^(x;yi,y2) = ln(ri + r2 + ria) , (B7) 

which satisfies 

^.9 = (B8) 
rir2 



in the sense of distribution theory. (We denote ri = |x — yi|, r2 = |x — y2|, ri2 = |yi — y2|-) 

Hj 



The essence of the present computation of /i^^'' is to consider a secondary kernel / defined 



by 



/(x; yi,y2) = ^(rrr2) 



ln(ri + r2 + r^) - ^ 



+ 7(^2^^! + ^2^2 - rir2) , (B9) 
o 



and whose main property is to satisfy 

AJ = 2g (BIO) 

in the sense of distributions. This property is easily checked on the expression ( [B9|) . We 
now recall that the potential Wij can be expressed in terms of the kernel g as 
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iy,,(x,t) = G' [[ rf3y^d3y^a(yi,t)a(y2,t)— — {(7(x;yi,y2)} . (Bll) 

Thus we see that a secondary potential defined by the same expression as (pil| ) but with g 
replaced by the secondary kernel /, i.e. 

w,ji^,t) =G^ ff rfVif^V2t^(yi,t)^(y2,t)— 7r7{/(x;yi,y2)} , (B12) 
necessarily satisfies 

/\w,j = 2Wij (B13) 
in the sense of distributions. Substituting ( pi3|) into ( p4D leads to 

= ^ STF,, FPb=o / (f^l^fU/^w.,, . (B14) 
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The next step is to integrate by parts the Laplace operator in the integrand of (|B14D . 
However, note that it is a priori not possible to perform this integration by parts ignoring 
the analytic continuation factor |x|-^ and the finite part symbol. For instance, the result 
(4.2) of Ref. [jlO[ invoked above does not apply to this case. We thus integrate by parts the 
integral (B14) keeping all terms coming from the differentiation of |x|^. This yields (using 
AU = -AnGa) 



riWl) 



— STFj,- / d y^awij + R 



(B15) 



where the first term has a compact support (we have removed on it the analytic continuation 
factor and the finite part symbol), and where the second term is explicitly given by 
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TiGc^ 
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B{B + 1) + 2Bx^dk 
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(B16) 



An equivalent expression for Rij which is convenient for our purpose is easily obtained by 
substitution of the expression ([B12|) of the potential Wij and use of t/ = / (i^y3cr(y3, t)/|x — 
ysl. We have 



STF,,. Ill rfVlC?V2rfV3^(yi)^(y2)a(2/3) ^^{^^yi'y2;y3)} (B17) 



where 

^(yi,y2;y3) = fPb=o 



5(5 - 1) - 2By, 



ci xlx 



B_2 /(x;yi,y2) 
|x- y3| 



(B18) 



To compute i^(yi, y2; y3), we need to control the pole part when i? ^ of the integral on 
the right-hand side of (pi8| ). This is because of the explicit factors B and B^ in front of 
the integral. [Note that it is important to keep the factors B and B^ in front since as we 
shall see the integral involves both a simple and a double pole when B — »• 0.] The pole part 
of the integral in ( [B18|) depends only on the behavior of the integrand at the upper bound 
|x| oo, so we need only to consider the asymptotic expansion of the kernel /(x;yi,y2) 
when |x| oo, or equivalently when yi,y2 — > 0. It is then easily shown that the only 
terms in the latter expansion of / which generate poles when are either of the type 

a regular solution of Laplace's equation, i.e. xl times a function of yi and y2, or of the 
type Xiln |x| times a function of yi, y2 (note that the expansion of / involves a logarithm 
of |x|). We shall slightly improperly refer to these terms as the "harmonic" terms in /. 
Their computation can be greatly simplified by noticing that in the asymptotic expansion of 
g = ln(ri -|-r2-|-ri2) when yi, y2 —>■ 0, only the three leading order contributions — constant, 
linear, and quadratic in yi, y2 — can contribute to /harmonic- Indeed, the higher order 
contributions, at least cubic in yi, y2, necessarily involve a function of x whose dimension 
is that of l/|x|" with n > 3 (because g is dimensionless) , and thus which will never yield 
a term of the type xl or x^ln |x| when multiplied by ri.r2 = x^ — x.yi — x.y2 + yi.y2- 
Computing the expansion of g when yi, y2 ^ 0, and then the expansion of /, we obtain 



/harmonic = ^yi .y2 ln(2 |x| ) - ^(x.yi 



x.y2j 



ln(2|x|) 



(B19) 
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[This computation of the "harmonic" terms in / bears a resemblance to the computation 
in Sect.IIIC of the "homogeneous" part of the function W.] The pole part of the integral 
in ( [B18| ) is then straightforwardly obtained from the replacements / — > /harmonic and |x — 
ysl"-*^ — > |x|^^(l + y3.x/|xp) (we do not need to expand |x — yal^^ further since /harmonic 
involves only terms with multipolarity or 1). We find a simple and a double pole. Then 
we apply to the latter pole part the operator B[B — 1 — 2y^d/dy^] present in ( P18|) (being 
careful about correctly handling the double pole), and compute the finite part at i? = 0. 
The result is 

^(yi,y2;y3) = ^yi.y2(ln2 + l) + i(yi.y3 + y2.y3)(ln2-^) . (B20) 
This is a non zero result; however the quantity of interest in ( pi7|) is 

dy\dy\ 



STFi,-^— ^{ir(yi,y2;y3)} = =^ i?., = , (B21) 



which shows that the contribution jfj^^^ given by (|B15| ) reduces to its first term, i.e. to the 
manifestly compact support integral 

= STF,, j d'xa w,, . (B22) 
This expression is our main result because the second contribution given by (|B5|) , is 



evaluated without problem by substituting into it the formulas ( PB] ) giving the divergences 
of the potential Wij and using at various places the function of Eq. ( |3.23| ) . Adding up all 
the terms constituting /|J^^'* to the first contribution liY^^ given by (|B22| ), we arrive finally 
at the following expression for the cubically nonlinear quadrupole term: 



= STF,, /// d^yid!'y,dh,a{yi)<y2Xy,) 



-Jvy r23 



dy\dy\ V3 



ln(ri3 + r23 + ryi) - ^ 
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^-r^{M - - y2yi) - + y\yi + yly'i) 

4ri3r23 6rf3ri2 



+ ^(%iy^ - Ay\y{ + hyy, - hy\yi) . (B23) 

0^13 J 

This expression is valid for a general fluid system. We can evaluate in an explicit way the 
differentiations with respect to yi and y2 it contains. In doing so one finds that thanks 
to the tracefree projection there is in fact no logarithmic term in ( [B23| ). After reduction 
of (|B23| ) to the 2-body case (using as usual (5-functions and formally discarding all infinite 
self-energy terms), we get 

= -^^p STF,, {(5mi + 2m2)2/i^' + (5m2 + 2im)y^ - Q{m, + m^)y\y{] . (B24) 
Changing notation with y\2 = y\ — y\ finally leads to 

= -^^p STF,,- {m,yY + m^y^ + 2(mi + m2)y^ + 2{rmy\ - m^yDyi^] , (B25) 
in complete agreement with the equation ( p.72| ) derived in the main text. 
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APPENDIX C: A COMPENDIUM OF FORMULAS FOR MOMENTS 



We list below the expressions of the time-derivatives of moments which are used in the 
computation of the waveform and energy loss. For the waveform: 
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For the energy loss: 
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The TT projections of relevant contractions of time-derivatives of moments with N as used 
in the waveform are: 
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- lb{Nnf{Nvfv}^ + b{Nnfv'^ - 4:{Nnfv'^ + 2,v^{Nnfn'^ 
^{Nvfn'^ -\v^n'^ + \v'^\\ . (C3i) 

Finally the "squares" of time-derivatives of moments as used in the energy loss are: 



(4f )2 = ?,2{vm) 



+ 



r 



1 - ^(212 -90i/) 



(130150 - 760071/ + 314421/^) ^ , 



2646 
26-^(1 -4t/) 



(C4a) 



15 

+ 480v^[30 - 7(241 - 122i/)] 

+ 18— -^^[420 - 7(4607 - 2074i/)] | 



126^^4^[21 - 27(146-611/)] 



/j(5) ^2 



512 

— (z/m) (1 - 3z/) — — 



450v^ + 578 



-v^ + 765 



Gm 



(/<?')^^(-)^^«^^ {1-^(17 -20.)}, 



128 GW 
— (i/m) (1 - 3i/) 



G^m^ 



2v^ + 2- 



Gm 



(C4b) 
(C4c) 
(C4d) 
(C4e) 
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